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Abstract 

In this paper we present a possible model of adaptive grids for nu- 
merical resolution of differential problems, using physical or geometri- 
cal properties, as viscosity or velocity gradient of a moving fluid. The 
relation between the values of grid step and these entities is based on 
the mathematical scheme offered by the model of scale- free networks, 
due to Barabasi, so that the step can be connected to the other vari- 
ables by a constitutive relation. Some examples and an application 
are discussed, showing that this approach can be further developed 
for treatment of more complex situations. 

1 A scale- free network model for adaptive 
grids 

We consider a differential problem of advective viscous flow in a domain Q, 
and let G a computational grid on Q with an initial grid step h. We consider 
a sistem of cartesian coordinates for the representation of the nodes. In ID 
the grid is a linear set of nodes, in 2D is a plane surface where nodes are 
vertices of squares with edge length equal to h, in 3D is a set of adjacent 
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cubes with squares of edge h as faces (see e.g. jTj, j7j). The grid is said 
adaptive if, during the evolution in time, its points are clustered in regions of 
high flow- field gradients, so that the step h can change (see e.g. [TJ). If x is 
the vector of the cartesian coordinates of a node and t is the time variable, 
then for an adaptive grid we have h = /i(x, t). 

We consider a real positive grid-field fconG defined on the nodes, i.e. if 
% is an index parameter of the i-th node and M is the total number of nodes, 
we have k = (ki)i<i<M = k(x.i,t)i<i<M- Assume that this field represents a 
physical or logical property associated to the flow, e.g. the number of fluid 
particles that interact by molecular forces with the particle on z'-th node, or 
the number of possible nodes where a particle can move from time t to time 
t + At using the physical rules of the flow. As other example, in the partic- 
ular matter of numerical weather prediction the topography of ground can 
influence heavly the computation of metheorological variables; its geometry 
is modelled by special grid using terrain-following coordinates (see e.g. [S]). 
Another interesting possible interpretation for ki is the number of values (u 3 -) 
of the velocity field, computed at time t on nodes of index j, that we could 
use in a numerical scheme for computing the value Ui on i-th node at time 
t + At with the desidered accuracy (see jS]). 

From these examples we can assume that k has high values in the regions 
of Q where the flow is turbulent or in general not laminar; in these cases, 
e.g., for a good accuracy of a numerical approximation of the velocity field 
u at time t + At, the number of values (uj), computed at time t, for com- 
puting the value Uj at i-th node can be greater then the number necessary 
for computing u at a node placed, at time t, in a region of laminar flow (see 
e.g. PP). Hence, where the grid step h has small values, the grid field k can 
have high values. 

We apply to field k the methods of a scale-free network, known in scientific 
literature as Barabdsi networks (0), assuming that the flow at time t + At 
on a node of the grid depends on the flow computed at time t at z'-th node 
with a probability function so defined: 



Using a continuum approach, Barabasi and Albert (|3|) assume that ki is 
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a continuous real variable and the rate at which it changes is proportional to 
n(fci): 



-of = ™n{h) (2) 

In |3] the factor of proportionality m is supposed in general constant. But 
if we want to describe an adaptive grid as a network changing in space and 
time, we suppose that the parameter m is in general a function depending 
on the node and on time variable: 



m = {mi)i<i< M = m(x i; t) (3) 

In this case the grid becomes a scale-free network of nodes with fitness, 
using the terminology of Bianconi and Barabasi ([3 ). Also, the dependence 
on the grid variable x suggests to consider the possibility of a law like (J2J) 
with the gradient V/cj at first member. 

The analytical dependence of m from Xj and t is determined by the physical 
and logical properties of the flow, such as viscosity, density, geometry of 
domain Q and the velocity field u itself. Hence in general we suppose that 
the mathematical expression of m is determined by a constitutive relation 
such as, e.g., 

m = m(u,Vu,/j,f) (4) 

where fi is the dynamical viscosity of the fluid. 

In the next section we'll present some possible examples of choices for the 
analytical shape of the function m(xj,t). 

Consider ki = /c(xj,t) = ki(t), m,i = m(xj,t) = nii(t). From (J2J) and the 
fact that by hypothesis the field k is positive, we can write 



dh m 



m , dt ( 5 ) 



If c = [c%^i))\<i<M is t ne se t of initial it = 0) values of k, we have 

rrii(s) 

and therefore 



log(ki) - log(d) = \ ^ M \ / ^ ds (6) 
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h = a exp ( / ™ l 8 ds J (7) 

The previous equation gives the relation between the value of the grid 
field k on the i-ih node and the variable parameter m. We can write it as a 
field equation in this way: 



k(xi,t) = c(xi) exp ( / ^ M " -ds ] (8) 



m(^j, s 

Now we state the relation between the variable grid step h = h(pc, t) and 
the grid field k. Let Xj and Xj two contiguous nodes of G, e.g. j = i + 1 in 
ID, or the nodes are two vertices of the same edge of a geometrical cube in 
3D. If hy = |xj — Xj| is the value of grid step between the nodes, according 
to the previous discussion we assume that hy depends on the inverse of the 
arithmetic mean of the grid field values: 

where A is a real positive function. If the grid is very fine, for every x G Q 
exists a topological neighbour ?7 X such that xi, x x G U x . In this case we can 
assume k t = kj and, using (JHJ), it is convenient to use the following definition: 

= Sfcij (10) 



2 Examples for adaptive grids 

In all the next applications we assume A constant and Ylj=i kj = & with S 
constant. We discuss some possible constitutive relation for the function m. 

Case 1. Let = for all nodes and for every t. Then, from (j7J), fol- 
lows that ki = Ci V£ > 0. Therefore, from ©, we have h^ = that is 
the grid is constant in time but it can be etherogeneous. In particular, if 
C{ = c > Vi, then the grid is homogeneous. This is the case appropriate for 
a laminar flow, for which, in general, an adaptive grid is not strictly necessary. 
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Case 2. Let m ; = m ^ 0. In this case 



mt 

h = Ci e s (11) 

and 



/iij = e"— (12) 

~t~ Cj 

Let m > 0. In this case the values of grid step at time t + At are smaller 
than the corresponding values at time t, so we can compute, as previously 
discussed, the velocity field of a flow which has a complexity growing in time. 
The flow is asintotically turbulent in all points of the domain Q. 
Let m < 0. In this case the values of grid step at time t + At are greater 
than the corresponding values at time t. The flow is asintotically laminar in 
all points of the domain Q. 

We can consider the particular case for which Cj = c > Vz and m << S. In 



this case e "s = 1 — therefore the grid step h is homogeneous in Q and 

S 



A ( mt , 
/> = - I 1 - I V/ e 





m 



(13) 



Case 3. According to adaptive mesh refinement methods (e.g. |lj), grid 
step values can strictly depend on flow velocity gradient Vu. Also, if the 
shape of domain Q is constant in time and the viscosity yu of the fluid is a 
variable in space (and in time), the grid G should be refined in the regions of 
turbulence, and in particular where /U has small values (0). For these reasons 
we can assume, as first approximation, that the parameter m is constant in 
time and depends only on the ratio under the hypothesis \x > 0. Using 
a formal Taylor expansion, we write 

/|Vu|\ |Vu| /|Vu|\ 2 
m = m + mi h m 2 (14) 

V fJ, J IX \ LL J 

where rrij are constants in all the grid G. 

If |Vu| << /i, then m = mo and we obtain the situations described in the 
previous cases. The flow becomes laminar or turbulent, depending on the 
sign of m . 

In the general case, assuming \i = /x(x) constant in time, from (jHJ) we have 

k = c exp ( — \— I \n 2 m + umilVul + m 2 |Vu| 2 l ds) (15) 
WS Jo 1 J 
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Let t > 0. Using the definitions ||Vu|| Ll ^ ^ = f Q |Vu(s)| ds, ||Vu||^ 2 ( 04 ) = 

J* | Vu(s)| 2 ds, and remembering that the coefficients nij are constants, we 
can expand to first order the exponential expression and obtain 

k = c + ^ + + (16) 

\ b fib fi 2 b I 

From ()10|) the grid step is 

h = — 7 2 V ( 17 ) 

c ( /i 2 ^ + /i 2 m t + fimi || Vu|| L1(0ji) + m 2 || Vu|| L2(0)t) J 

Consider a fixed t > 0. If Vu = 0, then the grid step for computing the 
solution at time t + At becomes 

h= J A - s (18) 

In general, if Vu ^ 0, for those regions in Q where at time t the viscosity /i 
has very small values, from (|T7j) follows that 

lim h = (19) 

This result is coherent with usual conditions of numerical stability for differ- 
ential problems where the advective terms are dominant (see e.g. [?]). 



3 Application to a diffusion-transport prob- 
lem 

In this section we show how, at least in a simple differential problem where 
the theorical solution is known, is possible to determine the grid field k. 

We consider the one-dimensional diffusion-transport problem ([7j) 

J -fxu"(x) + bu\x) =0 Vie (0, 1) . , 

{ u(0) = 0, u(l) = 1 { ' 

where \i and b are positive constants. If /i << b, case of dominant transport, 
the theorical solution is well approximated by the function 
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u(x) = e(^ (x " 1} ) (21) 
which presents a boundary layer effect at x = 1 (see Fig.l). 




Figure 1: Solution of the one-dimensional diffusion-transport problem in the case 
b = 100/x. 



For a good accuracy of a numerical resolution based on a possible numerical 
scheme, as Finite Differences or Finite Elements, we should impose a condi- 
tion on the local Peclet number Pe = that is Pe < 1. For the dominant 
transport situation, the condition implies that a constant grid step h could be 
very small, with large computational cost for the numerical resolution. From 
Fig.l we can see that the theorical solution has significant positive values 
only for x > 0.9, therefore for the resolution of the boundary layer effect we 
might construct a grid with a variable step h = h(x), where h is decreasing 
(0.8,1]. 



Because we want h ever smaller when x — > 1, on the basis of the Case 2 of 
the preceding section, we suppose that the entities m and ki are constant, 
m = mo an d ki = S> but from next considerations follows that should 
be sufficient the ratio be constant. Therefore, from (|7|). considering the 
unique variable x of the problem: 
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We suppose for simplicity q = c Vz, 1 < i < M. Then &i = k(0) = c, and 
from (jlUj) follows that h\ = — , so 

^) = — ferr ( 23 ) 

Let z > 1, the node from which we impose the condition Pe < 1. From 
definition of Peclet number we have 

bhi fmo \ 

< exp —x Vx > Xi (24) 

2fj, \ S J 

and, with no restriction, we can impose = exp (^a^) for the i-th node. 
In this way we obtain the required value for 



mo 
S 



Xi \ 2fJL J 



(25) 



For our aim, the argument of the logarithmic function must to be greater 
then 1, so the value of hi can be choosen so that hi> -r-. Note that, from 
(|2*n|) . - = so that we can consider the equation (|23|) as a constitutive 
relation of kind (J3J: 

VI = L log ( ^) (26) 



S Xi \ 2u' 

Using (|23j) . the adaptive grid nodes greater than be so determined 

(see Fig. 2): 

Xj + l = Xj H 7 7 r-r- (27) 

CTp (| iog (^)) 
In the numerical case of FigOU we have 

lim h(x) = 0.02, lim h = 0.0135 

Therefore the number of nodes for this heterogeneous grid is only few greater 
than the number of nodes for the homogeneous grid with h = 0.02 and, as 
shown in FigEl using the same numerical scheme the solution in the first 
case is more accurate than the solution in the latter. 
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Figure 2: A 2D representation of the adaptive grid for the one-dimensional 
diffusion-transport equation. The used values are hi = 0.1, b = 100^, Xi = 0.8. 



Conclusion and further developments 

We have developed a model of adaptive grid based on a constituive relation 
between the geometrical step and other physical variables of a differential 
problem. The basic mathematical hypothesis is a possibile relation con- 
necting the rate or gradient of the step and a probabilistic function which 
describes the physics or the geometry of the problem. 

Further developments can regard optimization of the technical scheme and 
theorical treatment of the accuracy of numerical solutions obtained using this 
type of adaptive grids. 
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Figure 3: Numerical solutions in the case b = 100/U, using Finite Differences 
method with decentrate approximation for the first derivative. 
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